Typeset with jpsj3.cls <ver.1 .0> 



Letter 



Fermi surface and mass enhancement in KFe 2 As 2 from de Haas- van Alphen effect measurements 
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We report on a band structure calculation and de Haas-van Alphen measurements of KFe2As2. Three 
cylindrical Fermi surfaces are found. Effective masses of electrons range from 6 to 18m,,, m e being the free 
electron mass. Remarkable discrepancies between the calculated and observed Fermi surface areas and 
the large mass enhancement (> 3) highlight the importance of electronic correlations in determining the 
electronic structures of iron pnicitide superconductors. 
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The discovery of superconductivity at T c — 26 K in 
LaFeAs(0, F) '' has given rise to intense experimental and the- 
oretical efforts to elucidate the superconducting pairing mech- 
anism and symmetry in iron pnictide superconductors (see 
Ref. 2 for a recent review). Since the development of real- 
istic theories of the mechanism requires detailed knowledge 
of the Fermi surface (FS), experimental determination of the 
FS is highly desirable. 

Accordingly, many angle-resolved photoemission spec- 
troscopy (ARPES) studies have been performed. 2 ' Their re- 
sults show some level of agreement in the FS and band disper- 
sion with conventional band structure calculations and mod- 
erate mass renormalization due to many-body effects. On the 
other hand, measurements of de Haas-van Alphen (dHvA) 
or other quantum oscillations, which are bulk probes and al- 
low accurate determination of the FS cross sections and ef- 
fective masses m", are rather limited. dHvA measurements 
performed on the FeP compounds LaFePO 3 ' 4 ' and SrFe2P2 5 ' 
have shown that band shifts of up to ~0. 1 eV are necessary to 
bring band structure calculations into agreement with experi- 
ments and that the enhancement of effective masses over band 
ones is about two. Since high r c 's are found only in FeAs 
compounds, dHvA studies of FeAs compounds are more de- 
sired. However, because of the structural/magnetic phase tran- 
sitions, measurements on the alkaline-earth 122 parent com- 
pounds AFe2As2 (A = Ca, Sr, and Ba) 6 ~ 8 ' have observed only 
small FS pockets, which makes it difficult to draw an over- 
all picture of the electronic structures of these compounds. 
Very recently, dHvA measurements have been performed on 
BaFe2(Asi_ v P x )2 for 0.41 < x < I. 9 ' As one goes from x = 
1 to 0.41, where T c ~ 25 K, the electron FS's shrink and the 
mass enhancement factor increases from ~2 to ~4. 

KFe2As2 is an end member of the high-r e binary alloy 
(Bai_ x Kj;)Fe2As2 with the ThCr2Si2 structure and has T c ~3 
K. 10, 1 The low-temperature resistivity exhibits a clear T 2 de- 
pendence with a large coefficient of A — 0.026 /iQcm/K 2 , 12 ' 
and specific heat measurements have found correspondingly 
large Sommerfeld coefficients: y exp = 69 or 93 mJ/K 2 mol- 
f.u. (f.u. = formula unit) for poly or single crystals, re- 



Table I. Experimental and calculated Fermi surface parameters. B || c. m e 
is the free electron mass. 
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spectively. 13 ' 14) These indicate the existence of moderately 
large electron correlations. 75 As nuclear quadrupole reso- 
nance measurements have shown that spin fluctuations (SF's) 
are much suppressed (compared with the optimally doped 
compound). 13 ' The first ARPES measurement 15 ' found a and 
P hole FS's at the Y point in the Brillouin zone (BZ) and e 
hole FS near the X point. The a and f3 FS's occupy 7 and 
22% of the BZ, while the e FS occupies 4 x 0.6% (because of 
the tetragonal symmetry, there are four e cylinders in the BZ). 
The authors assumed the double degeneracy of the j3 FS. The 
mass enhancement factor was found to be 2-4. A more re- 
cent measurement, 16 ' which was performed on single crystals 
from the same source as our dHvA samples, has confirmed 
the basic FS topology but with slightly different FS sizes, and 
has revealed that the a FS, not /?, is actually quasi-degenerate, 
i.e., two FS's with similar sizes exist. 

We here present results of a band structure calculation and 
dHvA measurements of KFe2As2. A substantial part of the 
large FS has been observed for the first time for the supercon- 
ducting FeAs compounds. Comparison with the band struc- 
ture calculation will reveal the importance of electronic cor- 
relations in the FeAs compounds. 

The electronic band structure of KFe2As2 was calculated 
within the local density approximation (LDA) by using a full 
potential LAPW (FLAPW) method. We used the program 
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Fig. 1. (color online), (a) Calculated electronic band structure, (b) density of states, and (c) Fermi surface of KFe2As2. Points of symmetry (solid circles) 
and lines of symmetry (open circles) in the Brillouin zone are explained in the top left figure of (c). 
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Fig. 3. (color online). Field-direction dependence of the Fourier transforms for the (a) (110) and (b) (010) planes, (c) is an enlarged view of (b) for a 
frequency range of a/j, . Note that the horizontal axis is F cos 8. A spectrum obtained at a field angle of 8 degrees is shifted vertically by 8 (in the units of 
the vertical axis). The amplitudes of some spectra are magnified (x5) or reduced (x0.5) for the clarity as indicated in the figure. The main spectra showing 
the frequency range < Fcosd < 5 kT are Fourier transforms over a field window 1 < B < 17.65 T, while the superimposed spectra for the range 2.9 
< F cos 8 < 4.5 kT are transforms over a narrower field window on the high field side, 10 < B < 17.65 T, so that the f peaks are enhanced. The dotted lines 
drawn for the a and e frequencies are based on Yamaji model, and the shading in (c) indicates a frequency spread due to a ± 1° distribution of the c axis 
orientation (see text). The lines drawn for the ( frequency are a guide to the eyes. No fundamental frequencies were observed above F cos 6 = 5 kT. 
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Fig. 2. (color online). Fourier transform of dHvA oscillations in inverse 
field at 0(oio> = 36.7°. The used field range is between 10 and 17.65 T. 
T = 0.08 K. dHvA frequencies are labeled with Greek letters. The peak 
marked "Cu(N)" is due to the Cu neck oscillation coming from Cu wire of 
the pick-up coil. A part of the raw oscillation data is shown in the inset. 
The oscillation data at high fields near 17 T may appear noisier than low 
fields, but it is not due to noise. It is because fast oscillations of the 3a/, 
frequency appear at high fields. 




Fig. 4. (color online). Effective masses versus magnetic field (B || c). Field 
ranges used for data points are indicated by horizontal bars with arrows at 
both ends. 



codes TSPACE I7) and KANSAI-06. The experimental crystal 
structure 18 ' including the atomic position zas of As was used 
for the calculation. Figure 1 shows the calculated band struc- 
ture, density of states (DOS), and FS. The calculated dHvA 
frequencies F and band masses intend for B || c are tabu- 
lated in Table I. The DOS at the Fermi level (£» is 58.4 
states/(Ry f.u.), which corresponds to the Sommerfeld coef- 
ficient of jband = 10.1 mJ/K 2 mol-f.u. As can be seen from 
the DOS plot, the states near Ep are mostly derived from the 
Fe 3d orbitals. Four bands 31-34 cross Ep, and the FS con- 
sists of three concentric hole cylinders at the F point of the 
BZ (bands-32, 33, and 34), small hole cylinders near the zone 
boundary (band-34), and a small hole pocket at Z (band-31). 
Although this FS is rather different from a previously calcu- 
lated one, 19 ' which has an electron cylinder at X instead of the 
hole cylinders around X, it is consistent with that determined 
by ARPES except that the small pocket at Z is not experimen- 
tally observed. 

Single crystals of KFe2As2 were grown from a potassium 
flux. Resistivity measurements on some of the grown crys- 
tals show the residual resistivity ratios of ~600. The dHvA 
measurements were performed in a dilution refrigerator and 
superconducting magnet by using the field modulation tech- 



nique. 20 ' The modulation frequency and amplitude are 67.1 
Hz and 10.4 mT, respectively. The field direction measured 
from the c axis is denoted by 8, and if necessary a subscript is 
attached to indicate the rotation plane. 

Figure 2 shows the Fourier transform of dHvA oscillations 
at f9(oio) = 36.7°. Six fundamental frequencies qj,, aij,, and 
are resolved. Figure 3 shows the field-direction dependence 
of the spectra. Note that the horizontal axis is F cos 8. Effec- 
tive masses of electrons were determined from the tempera- 
ture dependence of dHvA oscillation amplitudes between 0. 1 
and 0.43 K (0.2 K for &) for B \\ c (Table I). We also ex- 
amined the field dependence of the effective masses for the e 
and a frequencies; the masses are constant within the errors 
in the investigated field range as shown in Fig. 4. This is in 
stark contrast to a recent report of strong field dependence of 
the coefficient A of the resistivity. 21 ' 

In the case of a purely two-dimensional (2D) electronic 
structure, the FS would be a straight cylinder showing a single 
dHvA frequency F, and F cos 8 would be constant irrespec- 
tive of the field direction. In reality, there is a c-axis disper- 
sion of the electronic band energy, which leads to corrugation 
of the cylindrical FS and produces at least two dHvA frequen- 
cies corresponding to the maximal and minimal FS cross sec- 
tions. The simplest case where the c-axis dispersion takes a 
form of cos I c k z (I c is the interlayer distance) was considered 
by Yamaji, 22 ' and the angle dependence of the dHvA frequen- 
cies were derived. 

We assign each of the three pairs of the observed frequen- 
cies ei j,, a^h, and fy, to a corrugated cylindrical FS. The angle 
dependence of ey, is consistent with Yamaji model as shown 
by the dotted lines in Figs. 3(a) and (b). The angle depen- 
dence of a/j, can also be explained by Yamaji model if dis- 
tribution of the c-axis orientation in the sample is considered. 
The shading in Fig. 3(c) indicates the frequency spread cal- 
culated for the +1° distribution of the c-axis orientation. As 
the field is tilted from the c axis, the two frequencies or//, 
show a tendency to approach each other, but above around 
35° they split into small peaks, which spread out as expected 
from the calculated frequency spread. Despite the frequency 
spread, we notice an important characteristic of Yamaji model 
appearing at high angles; namely, the frequency peaks are en- 
hanced near 8 = +52.5°. This is because at this magic angle 
the maximum and minimum frequencies are expected to coin- 
cide and the whole region of the FS contributes to the single- 
frequency dHvA oscillation, resulting in an enhanced oscil- 
lation amplitude. The ( frequencies exhibit a distinctive an- 
gle dependence. The minimum frequency (i shows much less 
variation of F cos 8 than the maximum one This indicates 
that the corresponding FS cylinder is close to a straight one 
near the minimal cross section, and that it bulges out rather lo- 
cally around the maximal cross section. The effective masses 
of & and (i, differ considerably for B \\ c (Table I), but they 
become close as the field is tilted; m* = 11.1(4) and 12.5(7) 
m e , respectively, at %)io) = 36.7°. This also indicates that the 
modification of the FS is local. A similar but more extreme 
bulge has been observed for the larger hole FS of SrFe2P2- 5 ' 

In addition to the above frequency branches, many frequen- 
cies appear between 07 and (1 for field directions near B \\ c 
(Fig. 3). These frequencies can be expressed as F ai + nAF (n 
= 3, 4, 5, ...)> where AF = {F^ - F a/ )/l2. The frequencies 
n = 1 and 2 are not clearly resolved probably because they 
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are overlapped by the ay, frequency. Frequencies expressed as 
2F ai + nAF (n = 1, 2, 3, . . .) are also observed (small peaks 
appearing for F cos 8 > 4.6 kT in Fig. 3). The effective masses 
of the frequencies F a , + nAF range between those of ai and 
£i, and, as the field is decreased, their amplitudes are damped 
considerably faster than those of ai and £/. Therefore, they 
are most likely combination frequencies involving magnetic 
breakdowns. If the a and £ FS's are close enough for mag- 
netic breakdown to occur at some point in the BZ, there will 
be four such points because of the tetragonal symmetry. This 
can explain the frequencies with n = 3,6, and 9, and these 
frequencies are stronger than their respective neighboring fre- 
quencies. The origin of the observed three-times finer splitting 
however remains unclear. 

We identify our e, a, and ( cylinders with the e FS and the 
inner and the outer FS of the quasi-degenerate a FS observed 
by the ARPES measurements, 15 ' 16) respectively, and also with 
the band-34 small hole cylinder near X, the bands-32 and 33 
hole cylinders at T predicted by the band structure calcula- 
tion, respectively. The volumes of the observed cylinders are 
estimated from the averages of the minimum and maximum 
frequencies for B || c to be 1.1, 8.4, and 13% of the BZ, re- 
spectively. The total observed volume is (4x1. 1+8. 4+13=) 
26%. Since KFe2As2 is uncompensated, the total volume has 
to be 50%. We attribute the unobserved 24% to the B FS. 23) 
Using a 2D approximation, we can estimate the specific heat y 
from the measured effective masses. The observed FS's con- 
tribute 67mJ/K 2 mol-f.u to y (with a roughly estimated error 
of ±10%). This is to be compared with the measured single- 
crystal value of 93 mJ/K 2 mol-f.u., 14) and the difference is as- 
cribed to the B cylinder. 

Very recently, angle-dependent magnetoresistance oscilla- 
tions (AMRO's) have been observed in KFe 2 As 2 . 24) Two FS's 
are resolved and the estimated FS sizes are 12 and 17% of the 
BZ with a relative error of ±20%. They can be identified with 
the present a and f FS's. 

There are large discrepancies between the measured and 
calculated dHvA frequencies (Table I). As for the c-axis dis- 
persion, those of the a and f bands are considerably smaller 
and larger than calculated, respectively. In order to grasp 
how accurate band structure calculations usually are, we re- 
fer to our previous dHvA measurements on BaNi2?2 with the 
same ThCr2Si2 structure. 25 ' The band structure calculation for 
BaNi2P2 was done with the same procedure and codes as the 
present one. Discrepancies between the measured and calcu- 
lated dHvA frequencies for F > 0.5 kT were 7% at most and 
were typically much smaller. The remarkable discrepancies 
found in KFe2As2 suggest the necessity of treating electronic 
correlations beyond the LDA as suggested by previous the- 
oretical studies. Ref. 26 pointed out that self-energy effects 
arising from the spin-mediated interband interaction (mainly 
between electrons and holes) could cause band shifts rela- 
tive to LDA calculations, while Ref. 27 showed that when 
correlations were treated by the dynamical mean field the- 
ory the crystal field splitting of the Fe 3d orbitals were mod- 
ified, which would alter the FS. Since there is no electron FS 
and SF's are suppressed, 13 ' the latter seems more relevant to 
KFe 2 As 2 . 

The observed mass enhancements are about 3 to 7 for the 
a and f FS's. 28 ' They are broadly consistent with the specific- 
heat mass enhancement y eX plyband = 9.2. These enhance- 



ments are comparable to or even larger than the value of ~4 
found in the T c ~ 25 K compound BaFe2(Aso.59Po.4i)2- 9) Con- 
sidering the suppressed SF's and the field independence of the 
effective masses (Fig. 4), we do not ascribe the main origin of 
the mass enhancements specifically to low-energy SF's, but 
we ascribe it to band narrowing due to more general electronic 
correlations arising from the local Coulomb interaction on the 
Fe 3d shell. The importance of such high-energy correlations 
has been pointed out from recent optical studies on the iron 
pnictides. 29 ' It is also interesting to note a recent theoretical 
work, 30 ' which suggests that the mass enhancement increases 
with hole doping in iron pnictides. 

In summary, we have observed dHvA oscillations from the 
e hole cylinders near the X point of the BZ and the a and f 
hole cylinders at T. Effective masses of electrons range from 
6 to 18m e for B\\c, and the mass enhancements are 3 to 7 for 
the a and f cylinders. The specific-heat mass enhancement 
is 9.2. The unusual discrepancies between the observed and 
calculated FS areas and the large mass enhancements clearly 
demonstrate the necessity of including electronic correlations 
beyond the LDA in understanding the electronic structures of 
the iron pnictides superconductors. 
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